Energy loss of charged particles interacting with simple metal surfaces 
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Self-consistent calculations of the energy-loss spectra of charged particles moving near a plane- 
bounded free electron gas are reported. Energy-loss probabilities are obtained, within linear-response 
theory, from the knowledge of the density-response function of the inhomogeneous electron system. 
Self-consistent single-particle wave functions and energies are obtained by solving the Kohn-Sham 
equation of density-functional theory, and the electronic response is then computed either in the 
random-phase approximation or with the use of an adiabatic local-density approximation. Special 
emphasis is placed on the various contributions from collective and electron-hole excitations to the 
energy loss of charged particles moving parallel with the surface. The effect of the electronic selvage 
at a metal surface on the energy-loss spectra is also discussed, by comparing our full self-consistent 
calculations with those obtained for electron densities that drop abruptly to zero at the surface. 
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I. INTRODUCTION 

Charged particles interacting with metal surfaces cre- 
ate electron-hole pairs and, by virtue of the dynami- 
cally screened Coulomb interaction-bulk and surface col- 



lective excitations, i.e., plasmons. These excitations 
play a crucial role in the interpretation of surface elec- 
tron spectroscopy data, such as x-ray photoelectron spec- 
troscopy (XPS), Auger-electron spectroscopy (AES), and, 
reflection-electron energy-loss spectroscopy (REELS). El 
The interaction of swift electrons with surfaces has also 
attracted great interest in the field jpL scanning trans- 
mission electron microscopy (STEM)Jjl3 Equally, the in- 
teraction of moving ions with solids has -represented an 
active field of basic and applied physics and a great 
amount of research has recently been focused-pnJ;he case 
of ions that are incident at grazing angle. t3il3 Never- 
theless, existing calculations of energy-loss spectra in- 
voke either the local-dielectric, the hydrcjuiyiiamic, or the 
specular-reflexion model of the surface. L3 El An excep- 
tion is a recent self-consistent calculation of the stopping 
power of jellium jplanar surfaces for ions moving parallel 
with the surface.Ea 

In this paper, we present extensive self-consistent cal- 
culations of the energy-loss spectra of charged particles 
moving near a jellium surface. In the case of charged 
particles moving inside a solid, nonlinear effects are 
known to be c ruci al in the interpretation of energy-loss 
measurements;E3IHJ however, these corrections have been 
shown to be less impoctant when the charged particle 
moves outside the solid.Ej In Section II we present, within 
first-order perturbation [or, equivalently, linear-response] 
theory, general expressions for the energy-loss probabil- 
ity of charged particles moving along a definite trajec- 
tory in inhomogeneous media, and focus on the case of a 
bounded three-dimensional electron gas that is transla- 
tionally invariant in the plane of the surface. In Section 



III, we report the results of our full self-consistent calcu- 
lations of energy-loss spectra of charged particles moving 
parallel with the surface, which are found to satisfy sum- 
rules for particle-number conservation. Special emphasis 
is placed on the various contributions from collective and 
electron-hole excitations to the energy-loss probability. 
The effect of the electronic selvage at a metal surface on 
the energy-loss spectra is also discussed, by comparing 
our full self-consistent calculations with those obtained 
for electron densities that drop abruptly to zero at the 
surface. In Section IV our conclusions are presented. 

Unless otherwise is stated, we use atomic units 
throughout, i.e., = h = rUe = 1. 



II. THEORY 

We consider a recoiless particle of charge Zi mov- 
ing in an arbitrary inhomogeneous electron system at a 
given impact vector b with non-relativistic velocity v, 
for which [tetardation effects and radiation losses can be 
neglected.E3 Within first-order perturbation theory, the 
probability for the probe particle to transfer momentuHn 
q to the medium is given by the following expression :EZI 
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where L and A represent the normalization length and 
area, respectively, and VF(q, q';a;) is the screened inter- 
action 



VF(q,q';w) = J dr J dr' ^ 
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iy(r,r';cj) = u(r,r') + J dr^ J dr2 
xv{r, ri) x(ri, r2, cj) v{r2, r'). 
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Here, v{r,r') represents the bare Coulomb interaction 
and x(r, r',w) j^nthe so-called density-response function 
of the nicdium.E3 

Within a sclf-cnergy formalism, the decay rate of the 
probe particle is obtained from the knowledge of the 
imaginary part of the self-energy. In the GW approxi- 
mation, and replacing the probe-particle Green function, 
by that of a non-interacting recoiless particle, one finds:E£l 

= -2ZfJ2 JdrJ dr>*(r) 0}(r') 



xImM/(r, r', E, - Ej) 0,(r') (f>f{r), 



(4) 



where (/)i(r) represents the probe-particle initial state of 
energy Ei, and the sum is extended over a complete set 
of final states (f>f{r) of energy Ef. Describing the probe- 
particle initial and final states by plane waves in the di- 
rection of motion and a Dirac 5 function in the transverse 
direction, i.e.. 
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(5) 



where r_L represents the position vector perpendicular to 
the projectile velocity, one finds 



(6) 



T being the normalization time and Pq the probability 
for the probe particle to transfer momentum q to the 
medium, as obtained from Eq. (^. 

Alternatively, one may consider the energy that the 
probe particle looses per unit time due to electrpnic ex- 
citations in the medium. This can be written asE3 
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where p^^'(r, t) represents the probe-particle charge den- 
sity 



'(r,t) = Zi(5(r-b- vt). 



(8) 



and V {r,t) is the induced potential. To first order in 
p'^^^{r,t), i.e., within linear-response theory, one finds 
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x[W^(r,r',c.)-t;(r,r')]p^^*(r',i')- (9) 



Introduction of Eqs. (||) and (||) into Eq. (|^) yields the 
total energy lost by the particle 
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where Pq is, as in Eq. (|), the probability of Eq. (|l|) 
for the probe particle to transfer momentum q to the 
medium, and q • v represents the corresponding energy 
transfer. 

The results in Eqs. (|l|), (||), and (|l^) are general ex- 
pressions for the case of a classical trajectory in an ar- 
bitrary inhomogeneous electron system characterized by 
the density-response function x(r,r';w). In particular, 
in the case of a bounded three-dimensional electron gas 
that is translationally invariant in two directions, which 
we take to be normal to the z axis, the energy loss of Eq. 
( pO| ) may be expressed in terms of the two-dimensional 
Fourier transform of the screened interaction, as follows 
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where q|| and V|| are the momentum transfer and the ve- 
locity in the plane of the surface and z{t) represents the 
position of the projectile relative to the surface. Eq. ( pT| ) 
gives the energy that a charged particle moving with con- 
stant velocity along an arbitrary trajectory looses due to 
electronic excitations in an electron system that is trans- 
lationally invariant in two directions, as occurs in the 
case of a simple metal surface modeled by jellium. 



A. Parallel trajectory 

In the glancing incidence geometry ions penetrate into 
the material, they skim the outermost layer of the solid, 
and are then repelled by a repulsive ,|-acreened Coulomb 
potential, as discussed by GemmellH Through use of 
the appropriate effective potentials the ion trajectory z{t) 
can be calculated and the energy loss is then obtained 
from Eq. (|ll|). Here we restrict our attention to the 
case of charged particles moving with constant velocity 
V along a definite trajectory at a fixed distance z from a 
jellium surface, as approximately occurs under extreme 
grazing- incidence conditions. Eq. (nW then yields 
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where (—dE/dx) is the energy loss per unit path length 
of the projectile, i.e., the so-called stopping power of the 
electron system. 
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Eq. jl^ ) can be expressed in terms of P(g||,a;), which 
represents the probabiUty per unit time, unit wave num- 
ber and unit frequency for the probe particle to transfer 
momentum qn and energy to to the medium: 



dE 
dx 
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Alternatively, the stopping power of the system is of- 
ten described by means of P{iS)^ the total probability of 
exchanging energy lo with the medium: 
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where 
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dujLuP{Lu), (16) 



dq^l-mW{z,z;qii,u;), (17) 



with g|| = \/q'^ + (uj/v)^. 

The main ingredient in the evaluation of Eqs. ( p^ and 

( 



17|) is the screened interaction W{z, z'\ q\\ , a;). From Eq. 
, one easily finds 



W^(z,z';g||,cj) = w(z, z',9||) 4- j dzx J 

xv{z, zi; (j||) x(zi, Z2; t^) v{z2, z'; q\\), (18) 

where v{z,z';q\^) and x(z, z'; qn , w) are two-dimensional 
Fourier transforms of the bare Coulomb interaction and 
the density-response function, respectively. 

In particular, for z and z' coordinates that are well 
inside the solid, there is translational invariance in the 
direction normal to the surface and W {z , z' ; q\\ , to) can 
then be easily obtained as follows 

poo 

W{z,z';q^\,u;)^ dq. e'^^^^-^'Kiq) e-\q,iu), (19) 
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where q — \J q^i^ + q1 and e^^{q,Lo) represents the inverse 
dielectric function of a uniform electron gas, 



^{q,uj) = 1 -Ku((j)x((7,w), 



(20) 



v{q) and xi^.!^) being three-dimensional Fourier trans- 
forms of the bare Coulomb interaction and the density- 
response function, respectively. 

For z and z' coordinates that are far from the surface 
into the vacuum, where the electron density vanishes, Eq. 
(|l|) yields 



W{z,z']q\\,uj) = v{z,z';q\\) 



^e-«ii(-+^')5(g||,c.) 
911 



where g{q\\^uj) is the so-called surface-response function 

(7(g||,w) = -^ j dzi J dz2 e«li(^i+^^)x(^i,22;<Z||,^). 

(22) 

The|_energy-loss function l'nig{q\\,Lu) satisfies the / sum- 
ruleS 



dwwlm5(g||,w) = 2 7r^(7|| / e^'?"'^ n(z), (23) 



which applies to the case of a bounded three-dimensional 
electron gas whose exact density in the ground state is 
n{z). For z coordinates that are well inside the solid 
the electron density takes a constant value n, and for z 
coordinates that are far from the surface into the vacuum 
the electron density vanishes. 

In the long- wavelength limit {q\\ 0), 



(24) 



where lOs = ujp/\/2 and ujp — {AtttI) is the classical 
plasma frequency of a uniform electron gas of density n. 
Hence, in the q\\ — s- limit the energy loss is dominated 
by the excitation ofj^urface plasmons of energy LUg, as 
predicted by Ritchie .13 



B. The density-response function 

We consider a jellium slab of thickness a normal to the 
z axis, consisting of a fixed uniform positive background 
of density 



n+(z) 



n, —a < z < 
0, elsewhere. 



(25) 



plus a neutralizing cloud of interacting electrons of den- 
sity n(z). The positive-background charge density h is 
often expressed in terms of the Wigner radius r^, as 
l/n=(47r/3)r3. 

Time-dependent density-functional theory (TDDFT) 
shows that the exact density-response function of the 
electron system satisfies the integral equatior£3 



(21) 



X{z,z';qii,uj) ^ X°iz,z';q\\,uj) + j dzi j dz2 

xx°(z,z';g||,w) [v{zi,z2;q\\) + fxc{zi,z2;q\\,uj)] 

xx(22,z';g||,cj), (26) 

where x'^iz, z';q\^,uj) is the density-response function of 
non-interacting Kohn-Sham electrons 

X°{z,z';q\i,uj) = 2 ^(f>^{z)(f>*{z)(l)j{z')(p*{z') 
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and the kernel f^ciz, z';q^\,LL!) accounts for exchange- 
correlation (xc) effects beyond a time-dependent Hartree 
approximation. In Eq. (P7|), &{x) is the Heaviside step 
function, 77 is a positive infinitesimal, the energies Ei and 



Ej are 



and 



Ei — El 



(k|i+qii)' 
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(29) 



and the wave functions (j)i{z) and energies £i, which de- 
scribe motion normal to the surface, are the eigenfunc- 
tions and eigenvalues of the one-dimensional Kohn-Sham 
Hamiltonian 



1 cP 



(30) 



(p{z) being the electrostatic potential and Vxc{z) repre- 
senting the so^ealled xc potential of density-functional 
theory (DFT).EJ 

Within this scheme, the simplest possible approxima- 
tion is to neglect xc altogether and set the xc potential 
Vxc{z) and the xc kernel fxc{z, z'; q\^\ , uj) equal to zero. In 
this case, the one-dimensional single-particle wave func- 
tions and energies are the self-consistent eigenf unctions 
and eigenvalues of the one-electron Hartree Hamiltonian. 
The calculation of the density-response function is fur- 
ther simplified if the self-consistent electrostatic potential 
entering Eq. (BQ) is replaced by 



0, 



<y3(z) 



—a — zq < z < zq 



00, elsewhere. 



(31) 



where zq is chosen so as to ensure charge neutrality. This 
is the so-called infinite-barrier model (IBM).c3 Within 
this model, the one-electron wave functions are simply 
sines, and charge neutrality is easily found to yield 



zo = (3/16) Af [l + 0{\Fla)] 



(32) 



Exchange-correlation effects are usually introduced 
within the local-density approximation (LDA) of DFT, 
by replacing the xc potential at z by that of a uniform 
electron gas with the local density n{z). The xc kernel 
entering Eq. (|2^) is then set either cquaLto zero [this 
is the random-phase approximation (RPA)E3] or equal to 
the static {uj = 0) xc kernel 



dvxc{n) 



dn 



Siz^z'). (33) 



This is the sctcalled adiabatic local-density approxima- 
tion (ALDA).E3 

To compute the interacting density-response function 
x{z, z'] q\\ , Lu), we follow the method described in Rcf. 
pS. We first assume thatj-zt(z) vanishes at a distance 
zq from either jellium edge,Ej and expand the wave func- 
tions 4>i{z) in a Fourier sine series. We then introduce 
a double-cosine Fourier representation for the density- 
response function, and find explicit expressions for the 
screened interaction and the surface-response function in 
terms of the Fourier coefficients of the density-response 
function (see Appendix A). 

Great care was exercised to ensure that our slab calcu- 
lations are a faithful representation of the screened inter- 
action and the energy-loss probability in a semi-infinite 
medium. This issue is important, pin view of the sig- 
nificant quantum-size effects (QSE)c3 originated in the 
quantization of the energy levels normal to the surface: 
as the slab-thickness a increases new subbands for the z 
motion become occupied, thereby leading to oscillatory 
functions of a (the amplitude of these oscillations decays 
approximately linearly with a, and their period equals 
Xpf^, Xf = 27r/(37r^n)^/'^ being the Fermi wavelength). 
For each quantity a under study we considered three dif- 
ferent values of a. One such value is the threshold width 
a„ for which the nth subband for the z motion is first 
occupied. The other two values are a~ — a„ — Ai?/4 
and a+ = a„ + Ai?/4, and the infinite- width limit i s thi pn 
extrapolated with the use of the following relatioicilo 



a(a„ ) + a(a„) + a{a+) 



(34) 



n—n(z) 



Following this procedure to calculate the surface- 
response function g(q\\,uj), we have been able to prove 
that the sum-rule of Eq. (^) is satisfied for all values 
of q|| under consideration. The results presented below 
correspond to slabs with n — 12, for which a w 5 — 6Af- 



C. Simplified models 

For comparison, we also consider various simplified 
models for the screened interaction W{z, z';q\\,Lu) of a 
semi-infinite free-electron gas, which are all derived for 
electron densities that drop abruptly to zero at the sur- 
face. These are: (a) a classical model consisting of a semi- 
infinite medium of local dielectric function e(w), (b) 
classical and quantized hydrodynamic (HD) models^ 
and (c) a classical ipfimte barrier (GIB) or specular- 
reflexion (SR) modelp3'C3 which has the virtue of incor- 
porating dispersion effects by expressing the screened in- 
teraction in terms of the bulk dielectric function e{q,uj). 
Within these models and for z and z' coordinates that are 
outside the surface (z, z' > 0) the screened interaction is 
obtained through Eq. (^), from the knowledge of ap- 
proximate expressions for the surface-response function 



4 



1. Classical model 

Within tliis approach, the screened interaction is de- 
rived by imposing the ordinary boundary conditions of 
continuity of the potential and the normal component 
of the displacement vector at the surface [z = 0). For 
z,z' > 0, one then easily finds Eq. (|2l| ) with 



= 



e{uj) + l' 



(35) 



For a free-electron gas, the long-wavelength {q 0) 
dielectric function is 



eico) = 1 - 



U!(uj + irj) ' 



(36) 



and introduction of Eq. (^ into Eq. yields the 

long- wavelength limit of Eq. (|4 ) . Introducing this limit 
into either Eq. (^5|) or Eq. (17), one easily reproduces 
the classical expression of Echenique and Pendrylifl for 
the stopping power of a jellium surface, 



^^Zf'^Ko{2u;,z/v), 



dx ' 

where Kq is the zero-order modified Bessel function. For 



large values of z {z » v/lUs), Eq. (37) reduces to 



dE 



— = Z2 ^nuj./zv 



^ — 2 LJs z/v 



dx 



2v 



(38) 



to the energy-loss function lmg{q^^,uj) coming.^ from the 
excitation of either bulk or surface plasmons:C3 



(41) 



Img iqi\,uj) = - q\\ / dq^ S{lj-uj^ ) 



and 



Im5'^(9lh' 



2 9|| + 27^11 ujg 
and ijj^ 

of bulk and surface plasmons 



^Siu-coi), (42) 



respectively. Here, Loj^ and represent the dispersion 



(43) 



and 



{<)" 



-P^ql + Pqil y^2a;2+/32^2j _ (44) 



(3'') As in Eq. ^ q = .Iql + ql, and 



2/3 



(45) 



For the separate contributions to the sum-rule of Eq. 
( p3| ) coming from bulk and surface plasmons, integration 
of Eqs. dH) and (El) yields 



2. Hydrodynamic models 

In a HD model, the collective motion of electrons in an 
arbitrary inhomogeneous system is expressed in terms 
of the deviations from the equilibrium density. In a 
semiclassical approach, one writes and linearizes the ba- 
sic hydrodynamic equations, i.e., the continuity and the 
Bernuilli equation, and for a semi-infinite system finds 



iOZ. 



where 



A„ 



+ /32g2 _ (^(^ _|_ jj^) 



(39) 



(40) 



and (3 represents the speed of propagatioii-pf hydrody- 
namic disturbances in the electron system.Ell 

Within a quantized hydrodynamic model, one first 
linearizes the Hamiltonian of the hydrodynamic system 
with respect to the induced electron density, and then 
quantizes this Hamiltonian on the basis of the normal 
modes of oscillation, which are referred after quantiza- 
tion as bulk and surface plasmons. Hence, within this 
approach one can distinguish the separate contributions 



du uj Im^"^ {q\\,ijo) = 
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4 911 + 27, 
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and 



dujujlmg {q\i,uj)^- — uj 

4 g|| + 27,11 P 



(46) 



(47) 



respectively. It is then straightforward to show that for 
a semi-infinite system with a uniform electron density 
n that drops abruptly to zero at the surface the sum of 
bulk and surface contributions to the energy-loss function 
lmg{q\\,Lu) satisfies Eq. (p3|). 

In the limit as g|| — > the bulk contribution to the 
energy-loss function vanishes, and both the imaginary 
part of Eq. (^9[), on the one hand, and Eq. (|4^), on the 
other hand, yield the long- wavelength limit of Eq. (|2^). 



3. Specular reflexion model 

Either by neglecting, within the IBM, the intetference 
between incident and scattered electrons (CIBM),CJ or by 
simply assuming that electrons are specularly reflected at 
the surface (SRM)ll3, one finds 
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1 - es{q\\,uj) 
1 + es{q\\,uj)' 



where 



with q 



dqz 



(48) 



(49) 



+ and e ^{q,uj) being the inverse bulk 

dielectric function of Eq. (pO|). 

If dispersion effects are neglected altogether, thereby 
replacing the momentum-dependent dielectric function 
e{q,LLj) entering Eq. (|49| ) by a local dielectric function 
e{uj), Eq. ( ^ ) yields the classical prediction [Eq. (^]. 
Alternatively, if dispersion effects are incorporated in an 
approximated manner through the hydrodynamic dielec- 
tric function of a uniform electron gas, 



e(g,w) = 1 + 



(50) 



Eq. (^) is easily found to yie ld the hydrodynamic 
surface-response function [Eq. (p9[)]. 



III. RESULTS AND DISCUSSION 

We choose the bulk charge density n to be equal to 
the average electron density of valence electrons in alu- 
minum metal (r^ = 2.07), for which the Fermi momen- 
tum [qp = (Stt^tt,)^/'^] and bulk plasma frequency [up] are 
qp = 0.927 (oo is the Bohr radius, qq = 0.529 A) and 
iVp — 15.8 eV, respectively. We set Zi = ±1 and our 
results can then be used for arbitrary values of Zi , as the 
energy-loss probability is, within linear-response theory, 
proportional to Zf. 

In this section, we first show results for the energy-loss 
function lmW{z, z; q\^,uj) entering Eqs. ( p^ ) and (17). 
Fig. 1 shows lmW{z, z;qn,Lu), as a function of lu, with 



q|| = OAqp (for this small value of q\\ both bulk and 
surface plasmons are well defined excitations) and z co- 
ordinates that are either well inside the solid [z < —Xp] 
(Fig. la) or far from the surface into the vacuum [2 > Af] 
(Fig. lb). We have carried out self-consistent RPA and 
ALDA slab calculations for this quantity, as described in 
Sec. II. B, and have found the expected result that for 
these values of z they coincide with those obtained from 
Eqs. (Il) [z < -Xp] and (H) [2 > Xp]. Hence, the 
energy-loss function ImW{z, z;q\^,uj) represented in Fig. 
1 either does not depend on z (Fig. la) or depends on z 
through an overall factor of e"^"^!!^ (Fig- lb). 

In Fig. la, bulk-plasmon (dashed line) and electron- 
hole-pair (dotted line) contributions to the RPA energy- 
loss function lmW{z, z;q\^,Lu) are shown separately, as 
obtained from Eq. (p^, together with the total energy- 
loss function represented by a solid line. One sees that 
inside the solid the energy-loss spectrum is dominated, 
for small values of , by a continuum of bulk-plasmon 



excitations occurring at energies < cu < w^, where 
qc represents the critical momentum for which the bulk- 
plasmon dispersion enters the electron-hole-pair exci- 
tation spectrum. For Vs = 2.07 and = 0.4 qp, one finds 
= 17.6 eV and = 23.6 eV, and bulk plasmons can 
be excited by charged particles moving parallel with the 
surface with speed v > 1.13 wo (fo is the Bohr velocity, 
vo = 2.19 X lO^ms-i). 

For z coordinates that are outside the surface it 
had been generally believed that the energy loss orig- 
inates entirely in the iSxcitation of surface plasmons 
and electron- hole pairs .EJ Nevertheless, the continuum 
of bulk-plasmon excitations dominating the energy loss 
of charged particles moving inside the solid (see Fig. la) 
is still present for z coordinates outside the surface, as 
shown in Fig. lb, although the main contribution to the 
energy loss now clearly comes from the excitation of sur- 



face plasmons at uj 
one finds cj^^ 



9|| 

16.0 eVl. 



[for rs = 2.07 and q\\ = OAqp, 
The bulk-plasmon contribu- 
tion to Im W{z, z, g|| ,0;), as obtained for z = Xp within a 
quantized hydrodynamic model by introducing Eq. ( |4l| ) 
into Eq. ( |2l| ) (see also Ref. |4^), is represented in Fig. 
lb by a dotted line. The total SRM energy-loss func- 
tion, as obtained with the use of Eq. ( ^ ) and the RPA 
dielectric function e{q,uj), is represented in Fig. lb by 
a solid line. This curve shows that at low frequencies 
the energy-loss spectrum is dominated by the creation 
of electron-hole pairs, losses centered around cu^^^ are 
due to the excitation of surface plasmons exhibiting a 
finite linewidth, and bulk-plasmon excitations yield en- 
ergy losses at a; > that nearly coincide with the result 
one obtains within the quantized hydrodynamic model. 
We have also carried out self-consistent slab calculations 
of the energy-loss function liRg{q\\,uj), which for small 
values of (g|| < qp) has been found to be close to that 
obtained within the SRM. Nevertheless, differences have 
been found in the surface-plasmon energy , which 



shifts to lower frequencies as demonstrated earlierj£j 
also, the thickness of the slab is required to be very large 
for the self-consistent calculations to properly account for 
the high-energy spectrum originated in the excitation of 
bulk plasmons and electron-hole pairs. 

Fig. 2 shows self-consistent calculations of the prob- 
ability -P(w) for a charged particle to exchange energy 
u! with the medium, as obtained from Eq. (|l7|). The 
particle is assumed to move parallel to the surface with 
V — 2vo and two different z coordinates for which the 
electronic selvage at the surface is expected to play an 
important role: z — —Xp/2 (Fig. 2a) and z = Xp/2 (Fig. 
2b). Dashed, thick-solid and thin-solid lines give the re- 
sult of IBM, RPA and ALDA calculations, and dashed- 
dotted lines represent the SRM probability obtained with 
the RPA dielectric function e{q,uj). 

In the interior of the solid (Fig. 2a) , the bulk-plasmon 
energy-loss spectrum is known to be inhibited by the 
presence of surface-plasmon losses through the so-called 
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bregenzung or boundary effect predicted by Ritchie,i 
whose existence is due to the orthogonahty of the surface- 
plasmon and bulk-plasmon modes. Nevertheless, as the 
electronic selvage is changed from zero (SRM) to its 
actual structure (RPA and ALDA)j^ the creation of 
electron-hole pairs increases, the surface-plasmon peak 
diminishes, and a broad bulk-plasmon peak dominates 
the spectrum, showing that a proper treatment of the 
surface density-profile is crucial for the energy depen- 
dence of the energy-loss probability. These results are in 
qualitative agreement with the calculations reported in 
Ref. |2|, where the selvage structure is introduced in an 
approximate manner within the HD and the SR models. 

Bulk-plasmon losses occurring in the vacuum side of 
the surface (see Fig. lb) are not visible in the total 
energy-loss probability P(i^) which is dominated by the 
excitation of surface plasmons and electron- hole pairs, 
as shown in Fig. 2b. Also, this figure shows sub- 
stantial changes in the energy-loss probability as a re- 
alistic description of the surface response is considered, 
with an important shift of the surface-plasmon peak to- 
wards smaller energies, in agreement with the experi- 
mentalhi,|determined surface-plasmon energies of simple 
metals .EJO 

Figure 3 shows a surface plot of our full IBM (Fig. 
3a), RPA (Fig. 3b) and ALDA (Fig. 3c) calculations 
of the probability P{lo) for a charged particle moving 
with speed v ~ 2vo parallel to the surface. The plot is 
shown as a function of the energy loss w and the distance 
z from the particle trajectory to the surface. Although 
the energy-loss probability is found to be divided into 
losses centered around the bulk-plasmon energy (in the 
interior of the solid) and the surface-plasmon energy (out- 
side the solid) , this separation is not as clear as predicted 
with the use of simplified models for the surface response 
(see, e.g., Ref. |l]). As z — > — oo the energy-loss prob- 
ability P{uj) reaches a constant shape centered around 
the bulk-plasmon energy, which does not depend on the 
details of the electronic selvage at the surface and only 
depends on whether the xc kernel f^dz, z';q\\^uj) is set 
equal to zero (IBM and RPA) or not (ALDA). Outside 
the solid, the energy-loss probability, which is centered 
around the surface-plasmon energy, decreases with the 
distance z from the surface to the particle trajectory. 

Fig. 4a depicts our fuU IBM (dashed hue), RPA (thick- 
solid line) and ALDA (thin-solid line) calculations of the 
stopping power, as obtained from either Eq. (Q or ( |l^ ) 
as a function of z and with v — 2vq. In the interior of 
the solid, where the electron density is taken to be con- 
stant, both IBM and RPA stopping powers coincide with 
the well-known RPA stopping power of a uniform elec- 
tron gas. Short-range xc effects, included in the ALDA, 
provoke a reduction in the screening of electron-electron 
interactions, thereby increasing the energy loss. Outside 
the solid the electronic selvage at the surface plays a cru- 
cial role in the actual behaviour of the stopping power: a 
slow decrease of the electron density at the metal surface 
leads to a larger energy-loss probability [see also Fig. 2 



and Eq. (p3|)], and the IBM stopping power is, therefore, 
found to be too small. In the SRM the electron den- 
sity is assumed to drop abruptly to zero at the surface, 
which provokes a reduction in the electron-hole excitation 
probability, and the stopping power outside the solid (see 
Fig. 4b) is found to be even smaller than in the IBM. 
Low-energy excitations involve transitions from occupied 
electronic states near the Fermi level, which are sensitive 
to the actual density profile at the surface, and are found 
to play an important role in the energy-loss mechanism of 
charged particles moving with v < 2vq. Also plotted in 
Fig. 4b is the result of assuming that the stopping power 
for a charged particle that moves at a distance z from the 
surface can be approximated by thatj-of a uniform elec- 
tron gas with the local density n(2).Ej This often-used 
local-density approximation also yields an inaccurate de- 
scription of the position-dependent stopping power, due 
to the intrinsic nature of surface-induced excitations not 
present within this approach, and the results presented 
in Fig. 4 show the need for a self-consistent description 
of the surface response if one is to look at the energy loss 
of charged particles moving outside a solid surface. 

As the velocity increases the energy-loss spectrum of 
charged particles moving far from the surface into the 
vacuum is dominated by long-wavelength excitations and 
the stopping power is dictated by the integration of 
Lu Img((7|| , uj) , which as a result of particle conservation 
[see Eq. (|2^)] does not depend on the details of the ac- 
tual response of the solid. In this limit and with the aid 
of Eq. (|23) one easily finds the classical stopping power 
ofEq. 

The velocity dependence of the stopping power is 
shown in Fig. 5a, for a particle moving outside the sur- 
face at z = Xf- Our fuU IBM, RPA and ALDA calcula- 
tions are represented, as in Fig. 4a, by dashed, thick-sohd 
and thin-solid lines, respectively, and the SRM stopping 
power is represented by a dashed-dotted line, as obtained 
with the RPA dielectric function e(g,w). At low veloci- 
ties the energy-loss spectrum is dominated by intermedi- 
ate and short- wavelength excitations, even far from the 
surface into the vacuum, and a combination of the actual 
electronic selvage at the surface with the intrinsic nature 
of surface-induced excitations play an important role in 
increasing the energy-loss. At high velocities the energy- 
loss spectrum is dominated by the surface-plasmon exci- 



tation [see Eq. (24)] and all calculations converge with 
the classical limit of Eq. (37), as shown in Fig. 5b. As 
in Fig. 4b, the local-density approximation is also rep- 
resented in this figure, showing that this often-used ap- 
proximation cannot account for the energy loss originated 
in surface- induced excitations, not even at low velocities 
where the energy loss is entirely due to the excitation of 
electron- hole pairs. 
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IV. SUMMARY 



APPENDIX A: 



We have reported self-consistent calculations of the en- 
ergy loss spectra of charged particles moving parallel to 
a plane-bounded free-electron gas, in the framework of 
linear-response theory. 

We have found that the continuum of bulk-plasnion 
excitations dominating the energy loss of charged par- 
ticles moving inside the solid is still present for particle 
trajectories outside the surface. Nevertheless, these bulk- 
plasmon excitations are found not to be visible in the to- 
tal energy-loss probability -P(w) which outside the solid is 
clearly dominated by the excitation of surface plasmons 
and electron-hole pairs. 

As for the effect of the electronic selvage at the surface, 
we have found that the so-called bregenzung or bound- 
ary effect inside the solid is diminished, plasmon peaks 
are broadened, and the surface-plasmon peak is consid- 
erably shifted towards smaller energies. The electronic 
selvage at the surface has also been found to increase 
both the energy-loss probability and the stopping power 
for charged particles moving in the vacuum side of the 
surface. 

In the high-velocity limit and for charged particles 
moving far from the surface into the vacuum the actual 
stopping power is found to converge with the classical 
limit dictated by Eq. (p7|). However, at low and in- 
termediate velocities substantial changes in the stopping 
power have been observed as a realistic description of the 
surface response is considered, and we have concluded 
that a self-consistent description of the surface response 
is necessary if one is to look at the energy loss of charged 
particles moving outside a solid surface. Accurate mea- 
surements of the energy loss of protons being reflected 
from a variety jof salid surfaces at grazing incidence have 
been reported.Efl^E£l A theoretical description of these ex- 
periments requires that the ion trajectory z(t) be cal- 
culated and energy losses from the excitation of inner- 
shells be taken into account. Also, in real experiments 
band-structure effects might be important and the sur- 
face roughness might lead to additioiial energy loss due 
to the so-called Smith-Purcell effect .E3 Work in this di- 
rection is now in progress. 
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Here we give explicit expressions for the screened in- 
teraction and the surface-response function in terms of 
the Fourier coefficients of the density-response function. 

We first introduce the following double-cosine Fourier 
representation for the density-response function: 



oo oo 



X{z,z';q\\,uj) = ^ ^Xnin{q\\,(^) cos(m7rz) coa{mTz'), 

(Al) 



m— n— 



where d ^ a + 2zq a nd z — {z + zq + a)/d. 

Introducing Eq. (Al) into Eq. (18), we obtain the 
following expression for the screened interaction: 

OC CXD 

W{z, z';q\\,uj) = ^ ^ W"„„(q||,u;) cos(m7rz) cos(n7rz'). 



m— n— 



(A2) 



where 



~1p~ 



OC OC 



^ ^Vmm'{q\\)X7n'n'{q\\,U})Vn'n{q\\), (A3) 



M\\) 



m'— n' — 



qj + {rmr/dY 



2d 



Oil \l _ (^D^Q-Qwd] 

^ ^ (q2 + (n7r/d)2 



and 




(A4) 



(A5) 



Similarly, introduction of Eq. ( |A1[ ) into Eq. (|2 
yields: 

{<lh^): (A6) 



911 



m— n— 



where 



Mm^ll 1 - e^ll "^008(77171) 

d q?. + [imr/dY 



(A7) 
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FIG. 1. The energy-loss function, \mW{z,z;q\\,uj), as a 
function of u with = OAqp and z coordinates that are 
either well inside the solid [z < —\f\ or far from the sur- 
face into the vacuum [z > Xf]- (a) The solid line repre- 
sents the total RPA energy-loss function, as obtained from 
Eq. (Q with z < —Xf; dashed and dotted lines represent 
the corresponding bulk-plasmon and electron-hole-pair contri- 
butions, respectively, (b) The solid line represents the SRM 
energy-loss function, as obtained for z = Xf with the use of 
Eq. (Uq) and the RPA dielectric function €{q,u); the dotted 
line represents the bulk-plasmon contribution, as obtained by 
introducing Eq. into Eq. 

FIG. 2. The energy-loss probability P{io), as obtained from 
Eq. ( |r^ ) with V = 2vo and two different z coordinates: (a) 
z — — Af/2 and (b) z = Xf/2. Dashed-dotted, dashed, 
thick-solid and thin-solid lines represent SRM, IBM, RPA 
and ALDA calculations. The SRM probabilities have been 
obtained with the use of the RPA dielectric function e{q,ui). 
The damping parameter is taken to be 7 = ujp/10. 

FIG. 3. The energy-loss probability P{io), versus the en- 
ergy loss UJ and the z coordinate, as obtained for n = 2«o 
within the (a) IBM, (b) RPA and (c) ALDA. The sohd is in 
the region z < 0. 7 = cijp/lO. 



FIG. 4. Stopping power, as obtained from either Eq. ( |l4| ) 
or ( [1^ ) as a function of z and with v = 2vo. a) Dashed, 
thick-solid and thin-solid lines represent IBM, RPA and 
ALDA calculations, b) RPA (thick-solid hue), SRM -as ob- 
tained with the use of the RPA dielectric function e{q,ij)- 
(dotted line) and local-density (dashed line) calculations. The 
solid is in the region 2: < 0. 7 = ujp/10. 

FIG. 5. Stopping power, as obtained from either Eq. 
(^) or (|l^ ) as a function of v and with z = Xf/2. a) 
Dashed-dotted, dashed, thick-solid and thin-solid lines rep- 
resent SRM -as obtained with the use of the RPA dielec- 
tric function e{q, to)-, IBM, RPA and ALDA calculations, b) 
IBM (solid line) , SRM (dashed-dotted line) and local-density 
(dashed line). The classical prediction of Eq. (^) is repre- 
sented by a dotted line. 7 — ujp/100. 
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